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Abstract 

A reciprocating quantum refrigerator is studied with the purpose of determining the limitations 
of cooling to absolute zero. We find that if the energy spectrum of the working medium possesses 
an uncontrollable gap, then there is a minimum achievable temperature above zero. Such a gap, 
combined with a negligible amount of noise, prevents adiabatic following during the demagnetiza- 
tion stage which is the necessary condition for reaching T c — > 0. The refrigerator is based on an 
Otto cycle where the working medium is an interacting spin system with an energy gap. For this 
system the external control Hamiltonian does not commute with the internal interaction. As a 
result during the demagnetization and magnetization segments of the operating cycle the system 
cannot follow adiabatically the temporal change in the energy levels. We connect the nonadiabatic 
dynamics to quantum friction. An adiabatic measure is defined characterizing the rate of change 
of the Hamiltonian. Closed form solutions are found for a constant adiabatic measure for all the 
cycle segments. We have identified a family of quantized frictionless cycles with increasing cycle 
times. These cycles minimize the entropy production. Such frictionless cycles are able to cool 
to T c = 0. External noise on the controls eliminates these frictionless cycles. The influence of 
phase and amplitude noise on the demagnetization and magnetization segments is explicitly de- 
rived. An extensive numerical study of optimal cooling cycles was carried out which showed that 
at sufficiently low temperature the noise always dominated restricting the minimum temperature. 

PACS numbers: 05.70.Ln, 07.20.Pe 
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I. INTRODUCTION 



Reciprocating refrigerators operate by a working medium shuttling heat from the cold 
to the hot reservoir. The task is carried out by a controlled dynamical system. A change 
in the Hamiltonian of the system is accompanied by a change in the internal temperature. 
Upon contact with the cold side the temperature of the working medium is forced to be 
lower than T c -the cold bath temperature. A reciprocal relation is required on the hot side. 
Explicitly a quantum refrigerator is studied where the control of temperature is governed 
by manipulating the energy levels of the system. 

One of the main characterization of a refrigerator is the minimum temperature it can 
reach. The third law of thermodynamics already restricts this temperature to be the abso- 



lute zero 



a 



2]. Practically the minimum temperature is determined by the details of the 
mechanism of the heat pump. To investigate the cooling problem we study a model of a 
reciprocating quantum refrigerator. The present study is a comprehensive account following 
a brief version {3}. The main issues to be addressed are: 

• What are the restrictions imposed by the working medium? 

• What are the optimal conditions required to reach the minimum temperature? 

• Is there a minimum temperature above the absolute zero? 

To gain insight on these issues a reverse Otto cycle is considered where the working 
medium consists of interacting spin system. The magnetization/demagnetization stages are 
carried out by varying an external magnetic field which alters the energy levels of the work- 
ing medium. Such a model is a simplified version of adiabatic demagnetization refrigerator 
(ADR) 0, 0, 0]. These refrigerators have found use in cooling detectors to very low tem- 
peratures in space missions but also in an attempt to replace the existing technology in 
2 1,8 !_§!]• m addition t 



home appliances 



refrigerators 
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the purpose of unraveling the relation between quantum mechanics and thermodynamics. 



The present paper fol 
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ows a series of studies on a first principle four stroke quantum engines 



34 



35 



361 ] , where it was demonstrated that the model engines displays 



the irreversible characteristics of common engines operating in finite time 37]. 



2 



A generic working medium possesses a Hamiltonian that is only partially controlled ex- 
ternally: 

H = H int + K ext (u) (1) 

where ui = uj{t) is the time dependent external control field. Typically, [Hj nt , fl e xt] 7^ 0., 
therefore, [H(£),H(£')] ^ as a result a state diagonal in the temporary energy eigenstates 
cannot follow adiabatically the changes due to the control. The inability of the state to 



follow the change in the energy spectrum is source of quantum friction |32j, [33j, [3J, [35, 
361 ] . This friction limits the performance of the heat engine as well as the heat pump. 
There is an intimate connection between adiabatic following and the ability to reach cold 
temperatures. Since friction limits the performance almost perfect adiabaticity is the key 
to low temperature refrigeration. In this study we will explore the prospects of almost 
frictionless refrigeration cycles. Typically, the internal interaction in the working medium 
leads to an uncontrollable finite gap J in the energy level spectrum between the ground and 
first excited state. We will show that this gap combined with unavoidable quantum friction 
will be linked to a finite minimal temperature. 

A good characterization of the deviation from adiabticity is the difference between the 
von Neumann entropy of the state S vn = — tr{p log p} and the energy entropy defined by the 
projections on the energy eigenstate Se = — Yl Pj ^°SPji where pj is the population of energy 
state j. Equality is obtained only for perfect adiabatic following and thermal equilibrium 



33]. 



The present study explores the properties of quantum first principle four stroke heat 
pumps, with emphasis on the approach to absolute zero. In a previous study based on a 



phenomenological heat pump 31] we have found a linear relation between the cooling rate 
and the cold bath temperature T c . Does this relation survive when first principle quantum 
dynamical consideration are accounted for? 

II. THE CYCLE OF OPERATION, THE QUANTUM HEAT PUMP 

The working medium in the present study is composed of an interacting spin system. Eq. 
([I]) is modeled by the SU(2) algebra of operators. We can realize the model by a system of 



3 



two coupled spins where the internal interaction is described by: 

Hint = \hJ (&l ® & 2 X - ar\ ® & 2 y ) = HJB 2 (2) 

where <x represents the spin-Pauli operators, and J scales the strength of the inter particle 
interaction. For J — > 0, the system approaches a working medium with noninteracting atoms 



3l| . The external Hamiltonian represents interaction of spins with an external magnetic 
field: 

H-ext = -Ku{t) (&] ® i 2 + i 1 ® tr 2 ) = w(t)Bi . (3) 

The SU(2) is closed with B 3 = ± (<xj ® <r^ + o-* ® ct 2 , ) and [Bi, B 2 ] = 2zB 3 . 
The total Hamiltonian modeling Eq. (pQ) then becomes: 



H = n f w(t)B! + JB 2 J . (4) 
The adiabatic energy levels, the eigenvalues of H(£) are e\ = —hfl, e 2 /3 = 0, e± = hfl where 



Q = \/ijj 2 + J 2 . For J / there is a zero field splitting, an irreduceable gap between the 
ground and excited state levels. Eq. (jl]) contains the essential features of the Hamiltonian 
of magnetic materials 4J. 

The dynamics of the quantum thermodynamical observables are described by completely 



positive maps within the formulation of quantum open systems 38 



39 



40] . The dynamics 



is generated by the Liouville superoperator, C, studied in the Heisenberg picture, 

§ = i[H,A] +M A) + f • (5) 

where £^ is a generator of a completely positive Liouville super operator. 

The cycle studied is composed of two isomagnetic segments where the working medium 
is in contact with the cold/hot baths and the external control field uj is constant, termed 
isochores. In addition, there are two segments termed adiabats where the external field ui(t) 
varies and with it the energy level structure of the working medium. This cycle is a quantum 



analogue of the Otto cycle |32j]. Each segment is characterized by a quantum propagator 
U s . The propagator maps the initial state of the working medium to the final state on the 
relevant segment. The four strokes of the cycle in analogy with the Otto cycle (see Fig. []]) 
are: 

• hot isomagnetic (Isochore) A — ► B: the field is maintained constant uo = uo^ the 
working medium is in contact with the hot bath of temperature T^. leads to 



equilibrium with heat conductance T h , for a period of r h . The segment dynamics is 
described by the propagator U^. 

• demagnetization (expansion) adiabat B — > C: The field changes from ujh to u c in a 
time period of T\ IC . Cd = represents external noise in the controls. The propagator 
becomes Uh c which is the main subject of study. 

• cold isomagnetic (Isochore) C — > D: the field is maintained constant uj = uj c the 
working medium is in contact with the cold bath of temperature T c . C D leads to 
equilibrium with heat conductance T c , for a period of r c . The segment dynamics is 
described by the propagator U c . 

• Magnetization (compression) adiabat D — » A: The field changes from u c to Uh in a 
time period of T c h, £d = £n represents external noise in the controls. The propagator 
becomes U c h- 

The product of the four propagators, U s is the cycle propagator: 

Ucyc = UchUJAhcUh ■ (6) 

Eventually, independent of initial condition, after a few cycles, the working medium will 
reach a limit cycle characterized as an invariant eigenvector of U cyc with eigenvalue l(one) 



34| . The characteristics of the refrigerator are therefore extracted from the limit cycle. 



III. QUANTUM THERMODYNAMIC AL OBSERVABLES AND THEIR DY- 
NAMICS 

To facilitate the study of the dynamics of the cooling cycle we need a representation of 
the state p and the thermodynamical observables. The orthogonal set of time independent 
operators Bi, are closed to the dynamics. As a result they can supply a complete vector 
space to expand the propagators U and p. A thermodynamically oriented time dependent 
vector space which directly addresses the issue of adiabaticity is superior. This set includes 
the energy H and two additional orthogonal operators: 

H = w(t)Bj + JB 2 , L = -JBi + w(t)B 2 , C = fi(t)B 8 (7) 
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FIG. 1: Refrigerator cycle in the frequency entropy plane. The von Neumann entropy S vn = 
— tr{p\og p} (ABCD rectangle) as well as the energy entropy Se = — ^2p%logpi are shown [pi is 
the population of energy level i). The hot and cold isotherms are indicated. On the adiabats the 
energy level spacings change from to f2 c The demagnetization adiabat and the magnetization 
adiabats revolve exactly seven periods. On the isochores the energy level spacing remains constant 
and the entropy changes due to change in population. The cycle parameters are: J = 2, T c = 0.24, 
T h = 1.18, cj c = 0.1, uj h = 3.325, r c = 10.54, r h = 9.741, r hc = r ch = 12.81. 



To uniquely define the state of the system p the original set is supplemented with two 



f 



operators: V = Qf$4 = 
operator base the state p can be expanded as 



&l) and D 



cr\. With this 



P = 7* + 77 ( ( H ) H + ( L ) L + (C)C + (V)V + (D)D 



(8) 



4 ' n 

V and D commute with H. The equilibrium value of (V) is zero, and once it reaches 
equilibrium it does not change during the cycle dynamics. As a result the state p can be 
described in the energy representation by four expectation values: 

(l + ^D-E) U L + iC ) \ 




1 + ±(D + E)J 













1 - —D 



(9) 



where E — (H), L — (L), C — (C) and D = (D). From Eq. @ it is clear that when 
L = C = 0, p e is diagonal in the energy representation, then [p(t),H(t)] = 0. It therefore 
can be concluded that in complete adiabatic following L and C are maintained at zero value. 



A. The dynamics on the adiabats. 



In general the dynamics on the demagnetization adiabat is generated by C = Ch + £n 
where Ch = l[H, •] and H(£) is the time dependent Hamiltonian, Eq. (TJ1). The external 
noise generator C^ is defined later. The equations of motions for the dynamical observables 
Bi, B 2 , B 3 become: 
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■ > i 

d_ 
dt 
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B 2 
B 3 



(10) 



Our purpose is to evaluate the deviation from perfect adiabatic following. The equation of 
motion for the time dependent set H, L and C leading to the propagators Uh c and U c h are 
appropriate for this task. The integration to obtain Uhc and U c h will be carried out with 
respect to a new time variable dQ = Qdt: 

( TT \ /A -dm \ / H » 



d 
Qdt 



H 
L 
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Juj 










Juj 
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The ability of the working medium to follow the energy spectrum is defined by the adiabatic 
measure: 

We find that fi is a major parameter that characterizes the dynamics on the adiabats. When 
(jl = the propagator factorizes, the dynamics of H is independent of L and C. A large 
H will cause large non- adiabatic changes coupling H with L and C. We will show that 
constant /i minimizes the accumulated nonadiabatic transitions. 

Another nice feature of constant \i is that Eq. ffTTj) can be integrated leading to a closed 
form solution for the demagnetization and magnetization propagators Uhc and U c h- The 
consequence of stationary /i is a particular scheduling function of the external field u(t) 
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with time: 



Jf 



j 



(13) 



where f(t) is a linear function of time: fhc(t) = = (^ff- — j^J + j^. Swapping ft, 
for c in /(i) leads to the equivalent expression for the magnetization adiabat. 

The adiabatic parameter and the time allocated to the adiabat obey the reciprocal rela- 



tion: /j, hc 



The 



where K, 



he 



1 | Me 

j \ n c 



. Swapping c with h leads to fi c h and then 



K, 



eh 



-K, 



he ■ 



The solution is facilliated the time variable G, dO = Qdt. The final values of Qhc becomes: 



e 



he 



The 



he 



where: $ 



he 



arcsin(^) - arcsinf^M and > $ > — ~. 



n 



■ Q 



Eq. ( TTT]) is solved by noticing that the diagonal part is a unit matrix multiplied by a time 
dependent scalar. Therefore we seek a solution of the type Uh c = U1U2 where \U\,U%] = 0. 
The integral of the diagonal part of Eq. ( jTTj) becomes: 



(14) 



which can be interpreted as the scaling of the energy levels with the variation in Q. 

To integrate U 2 the non diagonal parts of Eq. (TTT1) . are diagonalized, leading to the 



eigenvalues 0, —i-Jq,iy/q, where q = \/l + /x 2 , and the propagator: 



/ l+/i 2 c us m(1~ c ) \ 



Uo 



\ 



/is 

1 

„2 



c _• 



(15) 



/ 



where s = sin(qQ) and c = cos(g6). The propagator U2 induces periodic mixing of H 
with L and C. As a result a diagonal p e Cf. Eq. (jHJ) will develop non diagonal terms. 
To characterize the deviation from perfect factorization of H from L and C, we define an 
adiabaticity measure 5 as: 

«J=l-Wf 1 W hc (l,l) . (16) 

where U^ 1 = ^-i is introduced to correct for the energy scaling. In the present context 
of noiseless dynamics and constant fi, 5 = 1—^2(1,1). When 5 = there is complete 
factorization. As will be described in Sec. [IV] 5^0 determines the minimum temperature. 
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The adiabatic limit is described by \i — > 0. Then Eq. (fTo"]) converges to the identity 
operator. These are the perfect adiabatic following conditions where 5 = 0. In general Eq. 
(115p describes a periodic motion of H L and C. Each period is defined by 

q O = 2ixl I = 0,1,2... (17) 

where I is the winding number. At the end of each period IA2 restores to the identity matrix. 
These are the periodic frictionless conditions where 5 = 0. For intermediate times (H) is 
always larger than the frictionless value 5 > 0. The amplitude of this periodic dynamics 
decreases when /i becomes smaller, Cf. £4(1,1) in Eq. ( ITST) . Constant /i is the minimum of 
5 Eq. (HED (Cf. appendix El 

The frictionless conditions define a quantization condition for the adiabatic parameter /1: 

<$>hc, 

Examining Eq. (1181) we find that there is no solution for I = 0. The first frictionless solution 
I > \ leads to a minimum demagnetization time: 

\/ Z7T *— ' 



/' = ( if) -1) ■ dS) 




2tt^ 



r hc (mtn) = K hc \j [ —J - 1 . (19) 

From Eq. (fT9l) we can interpret that the minimal frictionless demagnetization time scales 
as Th c (min) oc 4, since it has a weak dependence on uj c and uo^ . The special closed form 
solution can be employed in a piecewise fashion to analyze other scheduling functions uj (t). In 
general we expect similar quantization of the solutions. The main observation of this section 
is that we can find families of periodic frictionless solutions where the energy restores to its 
adiabatic value every period. For fi — > these solutions coalesce with the adiabatic following 
solutions. Table I summarizes some of the notations used. 



B. The influence of noise 



Any realistic refrigerator is subject to noise on the external controls. The main point 
of this paper is that even an infinitesimal amount of noise will eliminate the frictionless 
solutions. The sensitivity to noise results from the requirement of precise control of the 
scheduling of the external field u(t). To observe this effect requires a model of the noise 
induced by the external controls. 

9 



First we consider a piecewise process controlling the scheduling of u in time. At every 
time interval, us is updated to its new value. For such a procedure random errors are expected 
in the duration of these time intervals described by the Liouville operator Cn- We model 
these errors as a Gaussian delta correlated noise. This process is mathematically equivalent 
to a dephasing process on the demagnetization adiabat (35!] . This stochastic dynamics can 
be modeled by a Gaussian semigroup with the generator 40j, l4lj : 



|[H,[H,A] 



(20) 



d 

Tut 



L 



which is termed phase noise, Eq. ( 12"U]) . An equivalent dynamics to Eq. ([2"01 is also obtained 
in the limit of weak quantum measurement of the instantaneous energy [421 ] . For this noise 
model the modified equations of motion on the adiabats become: 

/ tt\ / n n \ / tt \ 

f 1-^ "I L • (21) 

The term 7 p describing the phase noise is assumed to be small. We therefore seek a pertur- 
bative solution: U a = lA^JU^U.-^ where the equations of motion U 3 can be obtained from the 



interaction representation: 



_d_ 

ttdt 



U 3 (t) = U 2 {-t) 



^ ^ 

--f p Q 
yO - lp Q J 



U 2 {t)U z {t) = W(t)U 3 (t) 



(22) 



where: 



W(t) = 7p fi(t) 



' £( a y* + 2(1 - c)) f^ 2 c+l) 
g(^ 2 c+l) ^ 



^(l-c)(/i 2 c+l)^ 

l-^(l-c) 2 ) 



(23) 



U3 describes the dynamics with respect to the reference provided by the unitary trajectory 
U.2- We seek an approximate solution for IA3 in the limit when fi —>■ 0, then U2 = I since this 
is the frictionless limit. Expanding Eq. ( 123]) to first order in /1 leads to: 



W(t) « 7p fi(t) 



11s /i(l — c) 

IMS 1 

\u(l-c) 1 



(24) 
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U-sijhc) is solved in two steps. First evaluating the propagator for one period of 0: for which 
Q(t) is almost constant, and then the global propagator becomes the product of the one 
period propagators for Z periods: ^3(17^) ~ ^3(6 = 2tt) 1 . The Magnus expansion {43] to 
second order is employed to obtain the one period propagator U^n): 



U 3 (Q = 2n) 



Mi+M 2 +- 



(25) 



where: M x = dQW(Q) and M 2 = \ f® dQdQ'[W(9), W(6')] + ...) The first order 
Magnus term leads to: 

/ 



W 3 (© = 2tt) 



Mi 



1 n{l - e ~ 2 ^ n ) \ 

e -2^7 P ^ 

\ - e' 27T ^ n ) e -^j P n 

which to first order in /x, 5 stays zero. U3 does not couple (H) with (L) and (C). 
The second order Magnus approximation leads to: 



(26) 



/ 



( c -s ^ 



W 3 (6 = 2rr] 



M 2 



SCO 







(27) 



V 



where S = sin a and C = cos a. a = 7 p f27r/i-y/9/i 2 + 4 and as /i 2 <C |, 



$ft c 7pfiy, Cf. Eq. ( JTBl . The condition fx 2 <C | can be transformed to / ^> An adiabat 
with a small number of revolutions Z ~ 10 already fulfills this condition. We now combine 
the second order propagator £4(t/j C ), for Z revolutions. It has also the structure of a rotation 
matrix identical to Eq. f[2"T|) . with a new rotation angle a = ai, where: 

r {Q h + uj h ){Q c -uj c 



2-K r ) p n I Vt{Q)dQ = 7T7 p Jln 
Jo 



(28) 



The asymptotic value of ati is finite when /i — > 0. 

For the quantization conditions when 1A% — 1 the deviation of W3 from the identity 
operator defines 5. Asymptotically as li — > and c<j c <C J, 



5min = 1 - COs(«i) « TT^^ln^/J] 



(29) 



Another source of external noise is induced by fluctuations in the frequency u(t). Such 
a term represent Markovian random fluctuations in the external magnetic filed. If the 
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fluctuations are fast compared to 2tt/Q, such noise can be described by the Lindblad term: 
C UJ X = - la u 2 [B 1 ,[B u X]]. 



( J 2 Jlu n \ 



_d_ 

ttdt 



U 3 {t) 



L) 



a 2 n 2 u 







Jul Ul 

n 2 n 2 
1 



U 2 {t)U 3 {t) = W(t)U 3 (t) 



(30) 



We seek an approximate solution for the quasistatic limit when \i — > 0. Expanding W in 
Eq. ([30l to zero order in /i leads to: 

\ 



^ J 2 Jujc —Jus 



W(0) 



-la 



n 3 



\ 



Juc uj 2 + s 2 J 2 scJ 2 
-Jus scJ 2 J 2 c 2 + u 2 



(31) 



/ 



We calculating the propagator for an integer number of periods the lowest order Magnus 
expansion becomes: W 3 (6 = 2vrZ) = exp (f* pU dQW(Q)^ , then the U 3 (l, 1) element 
decouples from the remaining part of the propagator and becomes: 



W 3 (l, 1) = exp 



Jo 



fi 3 (6) 



(32) 



Eq. (132!) can be integrated and since IA2 = 1 for an integer number of revolutions then: 



<5=1-Z4(1,1) ~ \- e la ^ Thc . (33) 



The smallest 5 is achieved for a one period cycle, Eq. ffTUl) then: 5 min w la J -qt ■ The phase 
noise and the amplitude noise have a reciprocal relation with respect to I. Phase noise is 
maximized for small I and amplitude noise for large I. Another possible source of noise is 
caused by fluctuation in the interaction energy H-mt- Analysis shows that such noise will 
lead to a similar expression to Eq. ( 1331) where J 2 is replaced by u 2 . 



C. The dynamics on the isomagnetic segments (isochores) 



On the isochores the equation of motion lead to equilibration with the hot and cold baths 
respectively. During the process the Hamiltonian is constant which leads to a factorization 
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TABLE I: Notation and definitions 
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(34) 



where T 



nn 

e fc s T and i£, 



eq 



MI(k + - k-)/T = hn(e k B? - 1)/Z, where 



Z = 1 + 2e fc s T + e 's T . In Eq. (El) H decouples from L and C. 
Eq. ( 134T) can be integrated leading to: 

e" r *(H(0) 



H(f) 
L(t) 
C(t) 



e" r * (L(0) cos fit - C(0) sin fit) 
e~ rt (C(0) cos fit + L(0) sin fit) 



(35) 



D(t) = D(0)e- 2r * + i(H(0)E eg (e- r '-e 
From Eq. (1351) the propagators U c and Uh can be constructed. 



E 2 eq (e~ rt 
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IV. THERMODYNAMIC AL RELATIONS 



The maximal efficiency r) max , of a heat engine is limited by the second law to the Carnot 
efficiency. For the quantum Otto type cycle the efficiency is limited by the ratio of the 



44j. As a result we obtain the series of 



energy level difference in the hot and cold sides 33 
inequalities: 

V max = 1 - ^ < 1 - — < 1 - ^ (36) 

' 'ft T h 

In the operation as a refrigerator the inequality in Eq. (|36|) is reversed. This imposes a 
restriction on the minimum cold bath temperature T c : 

T c > ^T h , (37) 
"ft 

fl c is limited by J and for the limit Uh > Jwe obtain: 

T c > —T h , (38) 

On the cold side the necessary condition for refrigeration is that the internal energy of the 
working medium at the end of the demagnetization is smaller than the equilibrium energy 
with the cold bath (Cf. Fig. [Hand ED- 

(H) c < (H) e ,(T c ) = -Ml c (l - 2e fc ^J , (39) 

where (H) e? (T c ) is approximated by the low temperature limit Ml c ^> ksT c . On the hot 

isochore the lowest energy point B, that can be obtained, is in equilibrium with T^. 

( h(l h \ 

(H)s > (H)eq(T h ) = -Wt h ( l-2e ^ J . Under these conditions L = C = 0. The 
change in (H) in the demagnetization adiabat leads to: 

(H)c « (1 _ 5) (H) B . (40) 
"ft 

where 5 the deviation from frictionless solutions is defined in Eq. (flB"]) . Then the maximum 
heat that can be extracted per cycle becomes: 

Q c (max) = (H) e9 (T c ) - (H) c « 2M1 C (V^ - - ~5 J (41) 



M2>, 



The condition for refrigeration is Q c (max) > 0. When 5 e fe f> T h the minimum temperature 
becomes the Carnot limit Eq. ( 1381) . For sufficiently large a^, positive Q c (max) > leads to 

h(l c 

5 <2e k b Tc imposing a stronger restriction on the minimal temperature: 

Tc - -fc fl iog(<y/2) (42) 
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Due to the logarithmic dependence on the noise S the minimum temperature scales linearly 
with the energy gap hJ. Eq. (|42[) relates the minimum temperature to 5 the adiabticity 
parameter. 



V. POWER OPTIMIZATION 

The cooling power V c is the amount of heat extracted Q c divided by the cycle time r. 
For the frictionless solutions on the adiabats the heat extracted is obtained by considering 
the balance of heat and work required to close the cycle 30 



E h E c (e Xc — l)(e Xh — 1) / wk nn r . \ 

Q c = nn c (^ - -^ r 1 _ ) ^ h j « 2hn c ( c -*ft - c -**) f(x c , Xh ) , (43) 

where: x c = Y c r c and Xh = Y h r h . Optimizing the cooling power becomes equivalent to 
optimizing F ^ Xc ' Xh "> where t cvc = r h + r hc + r c + r ch is the total cycle time. For frictionless 
solutions the minimum time on the adiabats Th c and r c h is described in Eq. ( fl9l) . The optimal 
partitioning of the time allocation between the hot and cold isochores is obtained when: 

r h (cosh(r c r c ) - 1) = T c (cosh(T h T h ) - 1) . (44) 

When T/j = T c the optimal time allocations on the isochores becomes = r c . 

The total time allocation r = r iso + r adi is partitioned to the time on the adiabats r adi 
which is limited by the adiabatic condition, and the time Tj SO allocated to the isochores. 

Optimi zing the time allocation on the isochores subject to ( |44l) leads to the optimal 



condition 36]: 



r c r C!/c (cosh(r /l r /l ) - 1) = smh(T h r h + T c t c ) - sinh(r c r c ) - smh(T h r h ) . (45) 
When T/j = T c = T this expression simplifies to: 

2x + Tr adi = 2 sinh(x) (46) 

(where x = T c t c = T h r h ). For small x Eq. ( 116]) can be solved leading to the optimal 

i 

time allocation on the isochores: t c = w (Tr a di/3) 3 /Y. Taking into consideration the 
restriction on the adiabatic condition this time can be estimated to be: r c = oc ^ (j) 3 - 

We can now expect two limits for the optimal cooling power the first when Y is sufficiently 
large the cycle time r cyc will be dominated by the time on the adiabats then for large uj c 



V c {max) oc hJ 2 e k t T - (47) 
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When the heat transfer time dominates, r c > 77^ then: 

V c {max) oc h—^e"^ (48) 

Noise on the adiabats modifies the optimal time allocation. Phase noise has its minimum 
for large values of I, Cf. Eq. (1291) . It approaches this minimum after a few revolutions 
independent of Kh c - The optimum power is a compromise between large time allocation on 
the adiabat to reach minimize noise and small cycle time to maximize power. As a result the 
scaling r adi oc 1/ J is still maintained, therefore Eq. fT4T|) or Eq. fj4*H|) will hold. For amplitude 
noise the minimum S is obtained for the minimum time frictionless solution which also leads 
to the scaling of power as in Eq. ( 1471) . 



VI. SIMULATING THE CYCLE 



After the segment propagators have been solved the cycle propagator can be assembled. 
For constant m the cycle propagator U cyc has a closed form solution. Other scheduling 
functions u(t) require numerical integration of the equation of motion Eq. fTTTT) . We have 
verified that our numerical integration coincides with the analytic expressions when available. 

The purpose of the simulation is to determine the optimal performance of the refrigerator. 
The cooling power was extracted from the limit cycle obtained by propagating the cycle 
iteratively from an initial state until convergence. The optimal cooling power was studied as 
a function of total cycle time r. For a fixed cycle time the heat extracted V c was optimized 
with respect to the time allocation on each segment. A random search procedure was used 
for this task. 

In general two types of cycles emerge classified according to the cycle time. The first 
are sudden cycles with very short periods r cyc <C ^ which are characterized by a global 
topology. These cycles are not presented and will be addressed separately [45(. The focus 
of the present study are cycles with a period comparable or longer than the internal time 
scale r cyc > 2tt/Q. The cycles of optimal cooling rate and minimum temperature are of this 
type. 

Fig. [5] and Fig. [3] present a typical cycle constructed with u{t) linear in t with optimal 
time allocation. Fig. [2] displays the entropy frequency plane. Fig. [3] shows the trajectory in 
the H, L and C coordinates. The positioning of the cycle with respect to the hot and cold 
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FIG. 2: Typical optimal cycle of refrigerator with linear scheduling, with T c > T™ n , in the 
1.(Q,Se), 2. (Q,Svn) plane (lower rectangle). The isotherms corresponding to the cold and hot 
baths temperatures, T c and are indicated. The difference between the energy entropy and the 
von Neuman entropy is the result of quantum friction: point A is higher than point D and point C 
is higher than point B. The parameters are: The cycle parameters are: J = 2, T c = 0.18, Th = 0.24, 
uj c = 0.1, uj h = 3.325, r c = 12.292, r h = 11.615, T hc = 18.016, T ch = 5.077. 

isotherms shows that it operates as a refrigerator with positive Q c . The end point D of the 
cold isochore is below the equilibrium point with the cold bath. On the scale of Fig. [2] this 
is hard to observe. One should also notice that Svn which is constant on the adiabats , and 
always a lower bound to Se, almost touches the minimal Se of the adiabats. The vertical 
distance from D to A, and from B to C is the result of quantum friction. 

The asymmetry between the demagnetization and magnetization adiabats can be noticed 
in both figure [2] and figure [3j The reason for this asymmetry is that the heat caused by 
friction in the magnetization adiabat can be dissipated to the hot bath. This is not true 
on the demagnetization adiabat where friction limits the possibility of heat extraction. This 
leads to very different time allocation > r^. The linear scheduling cycle should be 
compared to the cycle Fig. [T]and Fig H] where the friction is limited due to the quantization. 
The obvious difference is the symmetry between the demagnetization and magnetization 
adiabats. At the beginning and the end of the frictionless segments the von Neumann and 
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FIG. 3: Typical optimal cycle trajectory with linear scheduling shown in the H, L and C coordi- 
nates, for the same parameters as Fig. [2j Point A represents the beginning of the hot isochore. 
Point B represents the beginning of the demagnetization adiabat. Point C represents the beginning 
of the cold isochore. Point D represents the beginning of the magnetization adiabat. Notice the 
big difference between the demagnetization and magnetization adiabats. 

the energy entropy coincide. Periodic dynamics on the adiabats is also observed for optimal 
linear scheduling Cf. Fig El nevertheless frictionless solutions are not obtained. 

A. Numerical experiments 

We studied the optimal cooling cycles for a very large set of parameters for different 
scheduling functions. Fig [5] shows the optimal cooling power as a function of J/T c where the 
ratio 1Z = was maintained constant. The parameter 7Z addresses the distance of the 
operation conditions from the reversible limit where 1Z — 1. The simulations were performed 
for a predefined 7Z > 1 so that the second law is never violated Cf . Eq. ([551 . Fig. was 
obtained for a linear scheduling function of u(t) without the addition of noise. The sticking 
feature is that all graphs for different J values terminate at the same minimum J/T c . This 
graph has initiated our search for a possible explanation. In retrospect it represents the 
influence of uncontrolled numerical noise. Comparing to Eq. (1421) we can estimate the value 
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FIG. 4: Typical optimal cycle trajectory shown in the H, L and C coordinates, corresponding to 
Fig. Q] Notice that on this scale the isochores are barely observable. 

of 5 as ~ 10~ 7 . In addition all lines corresponding to different J values can be collapsed by 
shifting vertically by log J 2 . This finding shows consistency with the scaling of V c with J 2 
Cf. Eq. g7D. 

Simulations with constant m confirm the quantizations behavior of the optimal condi- 
tions. Figure [U] displays Q c for optimal cycles as a function of total cycle time r. The 
quantization of the cycle time is apparent corresponding to almost frictionless complete rev- 
olutions on the adiabats. The comb like function Q° p '(r) has a maximum at approximately 
I = 8 at the high temperature and I — 6 at the lowest temperature which is very close to 
the minimum temperature possible in these simulations. The maximum of Q" pt (r) is an 
indication of uncontrolled numerical noise in the simulation. Noiseless operation conditions 
would result in a flat comb distribution of Q° pt (r). 

Fig. [7J shows the entropy production as a function of cycle time. Only the cycles with 
very small entropy production corresponding to almost frictionless cycles, operate as refrig- 
erators. All non quantized cycles have a very large entropy production which decreases when 
the cycle time becomes longer. The quantized cycle (insert) have a very small entropy pro- 
duction which increases with the cycle time r. Fig. [7J demonstrates that quantum friction 
is accompanied by large entropy production. 
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FIG. 5: The logarithm of the optimal heat power V c as a function of J/T c for different J values 
(linear scheduling). On the graphs the ratio TZ = (T c f2^)/(T^O c ) = 1.453 is kept constant where 
oj c = 0.1 and T c /T h = 0.75. 

We attempted to identify the character of the numerical noise in the simulation. The 
procedure was to estimate the minimum temperature for a set of parameters J 1Z and 
C = Q c /Qh the magnetization ratio. Then we used Eq. (1411) to estimate 5. From the 
functional dependence of 5 on the parameters we tried to empirically asses the numerical 
noise in the simulation. In general we found both phase and amplitude noise. This can 
be observed in the trimming of both the high and low / ends of the comb in Fig. [61 In 
general we found that 5 decreased with C the compression ratio and with 1Z the deviation 
from reversibility. These dependencies were found for both constant m and linear scheduling 
where the constant m resulted consistently with a lower minimum temperature. The findings 
indicate that there is an additional source of numerical noise beyond the amplitude and phase 
noise. 

The existence of uncontrolled numerical noise hinders the study of the additional effects 
of the imposed phase and amplitude noise. The cycle simulation was repeated with the 
addition of phase noise Cf. Eq. (j2"Tj) . As can be seen in Fig. [H] an increasing amount of 
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FIG. 6: The optimal heat extracted Qc as a function of cycle time for three sets of temperatures: 
T c = 0.105, T h = 0.14 (Top: red) T c = 0.0975, T h = 0.13(Middle: purple) T c = 0.09, T h = 
0.12(Bottom: magenta). Results obtained by random search for stationary m with the restriction 
of r hc = T ch- Other parameters j = 2 uj c = 0.1, uJh = 3.32576 and no added noise. 

phase noise depresses Q c and moves the maximum to larger I, Cf. Eq. fl2Tj) . Numerical 
noise trims the high values of /. 

The quantization of the optimal cycles is independent of the specific scheduling. When 
the cold bath temperature T c is increased the quantization of Q c and V c is less pronounced. 
This can be seen in Fig |9] where the optimal power is plotted as a function of cycle time r for 
linear scheduling. The sharp comb structures of Fig. [6] is replaced by periodic modulation 
on top of a continuous background. At higher temperatures cycle with more friction can 
still operate as refrigerators and the quantization features are washed out. 



VII. SUMMARY 



Quantum friction is the result of the inability of the system to follow adiabatically the 
time dependent changes in the control Hamiltonian [H(t),H(t')] ^ 0. As a result the 
state of the system will develop non diagonal terms in the energy representation p e . The 
signature of this phenomena is an increase in the energy entropy Se- The key to cold 
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FIG. 7: The entropy production S u as a function of cycle time for three sets of temperatures 
corresponding to Fig. [6j Only the cycles with almost zero S u function as refrigerators. Their 
entropy production is shown in the insert on a very different scale. The quantum numbers are 
indicated. The break in the slope corresponds to the maximum V c in Fig. 



temperature refrigeration is no increase in energy in the demagnetization segment beyond 
the adiabatic limit. Perfect adiabatic following which requires infinite time will lead to 
frictionless demagnetization. Under conditions that the second law is fulfilled 1Z > 1 the 
cooling can continue to T c = 0. We then introduced an adiabatic measure \i to characterize 
the instantaneous nonadiabatic transition rate. The limit /i — > corresponds to perfect 
adiabatic following. Our first surprise was that constant fi led to closed form solutions for 
the dynamics. Moreover these solutions unraveled a quantized family of frictionless solutions 
5 = 0. These solutions are characterized by a state p commuting with the Hamiltonian at the 
beginning and end of the segment [H(0), p] = [H(t/j C ), p] = 0. Similar frictionless solutions 
were found for a working medium constructed from harmonic oscillators {46, l^j]. These 
frictionless solution can be carried out in a finite cycle time i.e. the cooling power does not 
vanish V c > 0. If these frictionless cycles could be realized they could operate to T c = 0. 

When we tried to simulate numerically the frictionless cycles we got into conflict. Any 
attempt resulted in a minimum temperature T c (min) which scaled linearly with the energy 
gap hJ. This observation eventually led us to the realization that any cycle is subject 
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FIG. 8: The optimal heat extracted Qc as a function of cycle time for three values of external 
phase noise 7 p : Top red 7 P = 0, middle orange j p = 10 -6 , middle magenta 7 P = 10 -5 , bottom blue 
7p = 2 • 10~ 5 . 7 P = 5 • 10 -5 did not result in positive Q c . Other parameters are as in Fig. [6l 

to noise. To follow this idea we constructed a model for external noise on the controls. 
Amplitude noise is the result of fluctuations in the magnitude of the external magnetic 
field. Since this noise term does not commute with the Hamiltonian it is not surprising 
that it will destroy the adiabticity, leading to 5 > 0. The surprise was the devastating 
effect of phase noise which commutes with the instantaneous Hamiltonian. Such a term 
can be the result of weak continuous measurement of energy during the adiabats. Such a 
measurement leads to partial collapse of the state to the energy representation. Naively 
one would expect this to lead to frictionless solutions. We have employed such an idea 
successfully to reduce friction in a quantum engine 35J. Nevertheless for a refrigeration 
cycle close to its minimum temperature phase noise accumulates leading to 5 > 0. Both 
types of noise are sufficient to eliminate frictionless solutions including the perfect infinite 
time adiabatic following frictionless cycle. 

Once the devastating effect of noise is appreciated it can be directly linked to a restriction 
on the minimum temperature. The minimum temperature T c mm depends on — 1/logJ Cf. 
Eq. (142p and will be on the order of the ener gy g ap hJ. This finding is consistent with 



experiments on demagnetization cooling of a gas 48| which obtained a minimum temperature 
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FIG. 9: The optimal cooling power Vc as a function of cycle time for linear scheduling. Upper 
plot (red) uj c = 0.5 and u>h = 3.9. Lower plot (blue) uj c = 0.1, Uh = 3.4. J = 1.25 T c = 0.725, 
T h = 0.966. 



an order of magnitude larger than the theoretical prediction [49| attributing the discrepancy 
to noise in the controls. Figure [10] shows the dependence of the minimum temperature of a 
refrigerator subject to phase and amplitude noise. The minimum temperature is related to 
the quantum number of the frictionless solutions I. The two types of noise show an opposite 
dependence on I. Amplitude noise favors small cycle times I = 1 while phase noise favors 
small m meaning large /. If both types of noise are operative the minimum temperature will 
be a compromise in Fig. [TU] at I = 25. Other sources of noise will also limit T c (min), for 
example our study was hindered by numerical noise. 

To conclude, it seems that any refrigerator constructed with a working medium possessing 
an uncontrolled energy gap will reach a minimum operating temperature on the order of the 
energy gap. For a working medium that has a controllable gap we found that if the gap is 
linear with T c there is no minimum temperature above T c = if the gap can be reduced to 
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FIG. 10: The minimum temperature as a function of the quantization number I. The diamonds 
represent phase noise and the circles amplitude noise. J = 2 
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APPENDIX A: OPTIMALITY OF CONSTANT jj, 

We show that constant /i is the minimum of the non adiabatic deviations i.e. minimum 
of of 5. We can transform Eq. ( TT21) to the differential equality: jidt = ^ leading to: 



The 



fi(t)dt 



duj 

Q3 



(Al) 



We decompose /i to a constant and a time dependent part /i = /U + ^ig(t). Without loss of 
generality we impose fio r hc = then J^ hc g(t)dt = 0. 

The first order correction to the propagator U2 due to time depdentce in /1 is the time 
average {l/r hc ) f^ hc dtU 2 (t). This will translate to a time average of 5. The dependence of 5 
Eq. (USD and Eq. OB on /i is: 



5 = fi 2 



1-c 



l + /i : 



2 ' 



(A2) 
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which is a monotonic increasing function of fi 2 with minimum at fi 2 = 0. The first order 
correction to U 2 will lead to 5 = 5 + 5i where S is the stationary result. Then expanding in 
fii will lead to: Si = (tq^) ^~ Sq H ° 9 2 (t)dt which is positive definite, therefore a stationary 
\x is a minimum of 5. 
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